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The like-sign dimuon charge asymmetry at the Tevatron: 
corrections from B meson fragmentation 
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The existing predictions for the like-sign dimuon charge asymmetry at the Tevatron are expressed 
in terms of parameters related to B mesons' mixing and inclusive production fractions. We show 
that in the realistic case when phase-space cuts are applied, the asymmetry depends also on the 
details of the production mechanism for the B mesons. In particular, it is sensitive to the difference 
in the fragmentation functions of B ( j and B® mesons. We estimate these fragmentation effects and 
find that they shift the theory prediction for this observable by approximately 10%. We also point 
out the approximately 20 % sensitivity of the asymmetry depending on which set of values for the 
B meson production fractions is used: as measured at the Z pole or at the Tevatron. The impact 
of these effects on the extraction of A| L from the D0 measurement is presented. 



I. INTRODUCTION 

Among the most natural places to look for deviations 
from the predictions of the Standard Model (SM) are CP- 
violation sensitive observables, since CP-violation in the 
SM is small. In a recent publication [![, the D0 collab- 
oration presented evidence for a 3.2a deviation from the 
SM. Together with the discrepancy in top-quark forward- 
backward asymmetry Q , this is the most significant cur- 
rent deviation from the SM. The implication of the mea- 
surement [l| is that the discrepancy might be due to a 
beyond the SM CP-violating effect in the neutral B me- 
son system. 

Given the intense interest in the results of Ref. [1| 
and their potential implication for beyond the SM (bSM) 
physics, in sec. IIIII we revisit the analysis with emphasis 
on B meson production effects characteristic for hadron 
colliders. To the best of our knowledge such effects have 
not been consider previously. In sec. IIII Dl we estimate 
the numerical impact of various higher order perturbativc 
and non-perturbative effects on the reported discrepancy. 

Accounting for S-production effects, in sec. II VI we con- 
sider the extraction of the flavor specific asymmetry ^4 SL . 
We observe that the inclusion of these effects further 
strengthens the expectation that both A| L and j4g L are 
deviating from their SM values. 



II. OVERVIEW OF THE MEASUREMENT 



The D0 collaboration measures [l[ the like-sign 
dimuon charge asymmetry A cxp , defined as 



_ N++ - N- 
A ™p = N+++N- 



(1) 



where N ++ (N ) represent the number of events con- 
taining at least two positively (resp. negatively) charged 
muons. A cxp contains both signal and backgrounds. 

The signal is defined as the same charge muon pairs 
originating from semileptonic decays of intermediate neu- 



tral B mesons: 



At 



(2) 



and is usually expressed as: 
b _ N(B°B° -> - N(B Q B° -> n~ yT X) 

As we will show in the following, the production mech- 
anism for the intermediate BB pairs does not decouple 
and also contributes to A b sl . 

In the context of Eq. the origin of the asymmetry 
is due to an oscillation B° <-> B° followed by a decay to 
a muon of the "wrong" sign Note that the inter- 

mediate neutral B mesons in Eq. ([2]) are a mixture of B® 
and B® mesons. 

An assumption built into the experimental analysis [l[ , 
and one we also adopt here, is that at parton level the 
production stage is symmetric in b and b. We will com- 
ment on it in the next section. 

The measurement involves the following cuts: < 2.2 
and px in the range 4.2GeV < px < 25GeV. Addition- 
ally, muons with px in the range 1.5GeV < pr < 4.2GeV 
are included if they have longitudinal momentum compo- 
nent \pz\ > 6.4GeV. The additional requirement \pz\ > 
6.4GeV for p T 's in the range 1.5GeV < p T < 4.2GeV is 
implemented by D0 in order to recover the acceptance 
for some of the muons in the forward region that would 
otherwise be lost due to the px cuts. 

There are two types of backgrounds contributing to 
Eq. (P). The first type is mostly due to Kaons (see Ref. [l[ 
for details); it contributes to Eq. Q]) but not to Eq. (|2|). 
The asymmetry in these events is due to a mechanism 
totally different from the one for the signal (essentially it 
is a detector effect). This type of background is roughly 
twice as large as the signal; it has to be subtracted from 
the data first. It is explained in |l| that this background 
is known with high precision. 

The second type of background consists of same sign 
muon pairs that contribute only to the denominators of 
Eqs. (HJ) and ((2]). The contribution from this type of 
background is relatively small [l[ ; see appendix [S] 
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III. THEORETICAL ANALYSIS 

The next step is to derive the SM prediction A t h for 
the quantity A b sl as defined in Eq. ([2]). In doing this we 
ignore the charge asymmetry A related to detection and 
identification of muons (see Ref. [l|) and assume it has 
been dealt with at the level of the data. We also ignore 
all backgrounds (discussed in appendix [A} . This way we 
consider only "signal" events that have un-equal proba- 
bilities for producing same-sign muon pairs of opposite 
charge. 

There are two such mechanisms. The first one is the 
signal described after Eq. ^fy where the same-sign muon 
pairs are resulting from the decay of intermediate B 
mesons after an oscillation. An example is the reaction 
b — > B — > /i~ ; b — >• B — >• B — >• . The asymmetry 
in this reaction is due to CP-violation in the neutral B 
meson system and is the subject of this study. 

As a second possibility, same-sign dimuon asymmetry 
can be generated from asymmetry in the production of 
bb pairs. A bb asymmetry cannot be generated pertur- 
batively within QCD, but at the non-perturbative level 
DGLAP evolution predicts non-zero asymmetry [6{ for 
the proton's bb parton distributions. While at the Teva- 
tron this should not produce any net asymmetry (since 
it is a pp machine) it might be of relevance for the LHC. 



A. Calculational framework 

The number of same-sign muon pairs N(pp — > B°B° — > 
H a H a X), with a = ±, can be derived within the frame- 
work of (heavy) di-hadron production at hadron colliders, 
including the decay of the produced hadrons. One has to 
calculate the partonic cross-section da(pp —> bb + X) dif- 
ferential in the rapidities and px of the observed quarks, 
as appropriate, and with full account of m\,. Then one 
has to convolute it with non-perturbative fragmentation 
functions describing the long-distance quark-meson tran- 
sition b — >• B q , q = d, s for each one of the two b, b quarks, 
followed by the decay B q — > \i + X. For more details 
about 6-production at hadron colliders see, for example, 
Refs. Adopting concise notations we have: 



iV ±=t (cuts) = / dPS da±± 

J cuts 



(3) 



with: 
da" 1 



da = 



E d °T 

i,j=d,s 



e d °r 

i,j=d,s 



D 



np _ , 






np 

b-+B s 






np _ ? 
6->flj 






np 

b^B-j ^ 







(4) 



(5) 



In Eq. ([3]) we have omitted an overall factor that cancels 
in Athj and <E>t (■ ■ • ) b denotes a convolution with respect 
to the variables of the parton b; similarly for b. 

The functions D np appearing in Eqs. (14151) are the non- 
perturbative fragmentation functions for the b — > B tran- 
sition that are detailed in sec. IIII Bl and D R ' W are the 
"right" and "wrong" B meson decay functions detailed 
in sec. UTTCl 

The functions da^ Tt in Eqs. ([11]) are the correspond- 
ing partonic cross-sections for producing a bb pair. Di- 
hadron production has been studied in the massless case 

una- 

It can also be applied to the massive case within 
the FONLL formalism of Ref. (HI, with the recent ad- 
vances of Ref. 14[ and with the good theoretical control 
over heavy flavor fragmentation |15l - [l8| . 

The intrinsic hard scale in the process is the invariant 
mass of the well-separated quark pair, which makes it 
possible to calculate observables even for small px 1 • 

As we have emphasized in the notation adopted in 
Eq. ([3]), the number of produced same-sign muon pairs 
depends on the cuts; the set of cuts is spelled out at the 
end of sec. HU Moreover, in presence of cuts, the effects 
of quark production, fragmentation and decay do not de- 
couple from each other - even after integration over the 
phase space. It is only in the limit of no cuts (i.e. un- 
restricted integration over the whole phase space) that 
such factorization might take place. 

Important assumptions have been built into 
Eqs. (|3I4I5[) . First, we have summed over all inter- 
mediate states which we take as B® and This is a 
standard assumption in works on the subject. Second, 
we have assumed that the production of the intermediate 
neutral B mesons factorizes from their decay. This is 
natural since the -Bd.s-mesons are well defined particles 
and different scales drive their production/decay. 



B. b — > B fragmentation functions 

The non-perturbative fragmentation functions D\^ B 
appearing in Eqs. (|4l5p are typically extracted from LEP 
data. We normalize them: 



D np - 

b^B a 



- f D np 

— Jqf, 



b^B a 



(6) 



where J* dxD^g — 1, and f q are measured production 
fractions. We introduce the following simplified notation: 



= D np 

b^B a 



= D f 



(7) 



The fractions f Sj d are measured at the Z pole and at 
the Tevatron. The two sets of measurements differ from 
each other by roughly 20%, although the fractions are 



1 Note that since the b-quarks are massive, mj, >> Aqcd> the 
hard scale is bounded from below by mj,. 
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generally assumed to be universal. The most up-to-date 
values can be found in Table 4 of Ref. [19| . Indeed, as long 
as the factorization theorem applies, they have to be pro- 
cess independent since they refer to the non-perturbative 
transition b — > B which takes place at scale ~ « Q. 
In the analysis of Ref. |l| the values measured at the 
Tevatron are used. 

The fragmentation functions Db s and Ds d are not 
known. The most accurately measured fragmenta- 
tion functions at the Z pole are for a combination of 
B ± ,B%, B° s mesons and A b [2(^22] ■ We expect that these 
two functions are different. First, the two fragmentation 
fractions are rather different (by approximately a factor 
of 4) which implies different non-perturbative dynamics 
in the two cases. Second, we can roughly estimate these 
functions by utilizing the fragmentation function of Pe- 
terson et al. [23| : 

D B ,(x)=Af-(l----^-) 2 ■ ( g ) 

X \ x 1 — X J 

where Af is fixed by the normalization and we take the 
parameter e q = m q /mf. In the following we will use 
rn b = 4.75GeV (as in Ref. Q) and m d = 70MeV,m s = 
170MeV. Our motivation for choosing these values for 
md, a is the following. The best value for e extracted 
in Ref. Q is e = 0.002, and it tends to decrease with 
the inclusion of higher order effects (see also 24] ) . This 
parameter is for the above mentioned combination of b- 
flavored hadrons. Assuming that the fragmentation func- 
tions for the charged mesons and baryons are softer than 
these for the B d s , our choice of mj jS leads to the follow- 
ing reasonable values for the fragmentation parameters 
e s = 0.0013, e d — 0.0002. A second constraint we have 
imposed on m^ s is that the difference m s — md is con- 
sistent with the mass difference of the corresponding B 
mesons (or of the current quark masses). 

In our subsequent analysis we need the fifth moments 
A = 5 of these functions — J* x N ~ 1 D Bq (x)dx for 

which we obtain D^] = 0.69 and D%\ = 0.82. In fact, 
the only information about the functions D B that will 
enter our analysis is their ratio 

k = D%]/D%1 =0.84. (9) 

We stress once again that our method of estimation 
should be adequate for its restricted use. Furthermore 
our estimates are consistent with the expectation for 
the difference between these fragmentation function dis- 
cussed in Ref. [251 ] . 

C. B o B oscillations and B — >• fi + X decay spectra 

There are two types of decays for the produced B 
mesons: "right" (R) decays B — > fi + + X (or B — > 
fi~ + X) and "wrong" (W) ones B — > (i~ + X (or 
B —> fi + + X). For semileptonic decays the "wrong" 



decays are due to an oscillation B o B followed by a 
"right" decay (see Refs. 5] for more information). The 
"right" decay is assumed to contain no direct CP vio- 
lation Q. Also, oscillation and decay are assumed to 
decouple from each other. Thus, the decay functions in 
Eqs. (jUS]) read: 

D%^ 1+ = D%^- = T(B q -+ B q )T q Dq , 

Dj q ^-=T(B q ^B q )T q Df, 

Dj q ^ + = T{B q B q )T q D* . (10) 

The decay function Dq in Eq. (|10p is normalized 

Jq 1 D^{x)dx — 1, and T q is the semileptonic decay width 
of B q . We have also made use of the following result 
about the time-integrated probabilities T(B q — > B q ) — 
T(B q B q ) [3]. The semileptonic widths for B° and B° 
mesons are almost equal T s « T^. In the following we 
will consider that the decay functions are also equal, i.e. 
= D d . We will not need their explicit form. 

D. Putting everything together 

Combining the results above, we get the following re- 
sult for the number of same-sign muon pairs: 

Aleuts) = ^fW. ( n ) 

i,j=s,d 

N+ 3 + = /,/';! M , / -»• Bi)T(Bj -> /i ; :/,,(cul>; , 
= ,/ :/>', -> Bi)T(Bj -> B,)^(cuts) . 

In the above equations we have introduced the follow- 
ing function: 

/W(cuts) = / dPSda% Ii ® b (DB i ®Df i ) b 

J cuts 

®- h (D Bj ®Df)- b ,(12) 

where /y(cuts) = /jj(cuts) and i,j — s,d. 
Finally, the result for A t h can be written as: 

N++ (cuts) -N— (cuts) = A num 
th N++(cuts) + N— (cuts) ~ A dcn ' [ ' 

where the function A num reads: 

A num = / s 2 r2/ ss (cuts)T(B s -> B S )T- 

+ f d r 2 d f dd (cuts)T(B d B d )T- 

+ f s fd^sTdfsd(cnts) 

x [T(B d B d )T- + T(B S B s )T d ] , (14) 

-4don = A uun i(T~->r+) , q = s,d, (15) 

and 

T± = T(B q -> B q ) ± T(B q B q ) . (16) 



4 



Eqs. (|13I14I15P are our main result. In general, the 
three functions f ss , fd d and f s d are different from each 
other and therefore Eq. (fl4|) does not factorize. Clearly, 
the function A t h becomes dependent on the experimental 
cuts, too. In order to gain more inside into the properties 
of this result, we consider limiting cases first. In the two 
cases: 

1. Db b = DB d and = for any cuts, 

2. fully inclusive case (i.e. integration over the full 
phase space) for arbitrary fragmentation and decay 
functions, 

we have f ss = fdd = fsd- In these two cases the functions 
fij factor out completely from Eqs. (|14I15[) and cancel in 
Eq. (fT3"|). Therefore, A t ^ becomes independent of the 
kinematics (i.e. cuts) and takes the usual form Q: 



ith 



f s T- + f d T- 
f s T+ + f d T+ 



in cases 1,2. 



(17) 



To get an insight into the general result Eqs. (|13|14|15|) . 
we consider the following simplification: 



/ij(cuts) w <f>i (cuts) 4>j (cuts) 



(18) 



i.e. the dependence on the two fragmentation channels 
factorizes. The motivation for this approximation follows 
the results of Ref. [7J : in the case of single B meson pro- 
duction in the same px range, the partonic cross-section 
behaves as a power da ~ 1 /pj, . In such a case the convo- 
lution becomes sensitive only to the n-th moment of the 
corresponding fragmentation function and the convolu- 



D 



(n) 



In the 



tion becomes a simple product, i.e. 4 
following we take n = 5. 

Within the approximation (|18p . the same sign dimuon 
asymmetry reads: 



Ah = 



/ s s (cuts)?7 + / d d (cuts)T d - 
/ s ^ s (cuts)T+ + / d ^(cuts)T+ 



(19) 



The above equation has the virtue that it exhibits the 
essence of the general result (|13I14I15[) (i.e. explicit de- 
pendence on cuts and production stage effects) but in a 
simplified manner that allows direct interpretation. In 
particular, we see that all the dependence on cuts and 
B meson production-stage effects can be absorbed into 
effective production fractions /° ff = /^(cuts). It is for 
this reason we believe that our analysis underscores the 
need for a better understanding of the difference in the 
production fractions /j measured at the Tevatron and at 
the Z pole. 2 

We next use the independence of (fT3"|) with respect 
to the normalization of the functions </>,• and divide the 



whole result by <pd 
sorbed in the ratio 



This way all </>-dependence is ab- 



D 



D 



(n) 



= K < 1 , 



(20) 



as follows from the discussion in sec. MI B[ and our de- 
fault choice is n = 5. 

Expressing the time integrated amplitudes T through 
the semileptonic asymmetries Ag h Q: T+ = Z q /T q , 
T- = A q SL Z q /T q , with Z q m 1 - 1/(1 + x 2 q ), we get: 

A t h = C S A^ L + CdAg L , 

Kf S Z s 



Kf S Z s + f d Z d 
fdZd 

nf s Z s + f d Z d 



(21) 



For consistency, the numerical values for the parame- 
ters / g (Tev) and x q are as in Ref. [l|, while fi(Z pole) 
are taken from Ref. [l9| . We stress again that within the 
approximation (|18[) all dependence on cuts and B meson 
production effects is contained in the parameter k. 



TABLE I: Ath as a function of the input parameters fi (Teva- 
tron versus Z pole) and k. Only uncertainties from fragmen- 
tation are shown. 



K 


1.00 


1.00 


0.84 


0.84 


fi measured at: 


Tev 


Z pole 


Tev 


Z pole 


A th x 10" 4 


-2.33 


-2.77 


Z - 04 -0.13 




(A th -^)/< [%) 
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9± 6 4 


28±1 



In table U we present A t h as a function of the param- 
eters fi and k. We present both its absolute value and 
the relative change with respect to the central value of 
4h = (-2-33t°-j)l(T 4 given in Ref. Q. We have stud- 
ied the sensitivity of our result due to variation of the 
fragmentation function parameters (i.e. md. s ) and the 
moment n of the fragmentation function. This is an ap- 
proximately ±5% effect, as indicated in table HI 

For comparison, the measured value [1] is yl cxp = 
(—96 ± 29)1CP 4 , where for simplicity, we have added the 
statistical and systematic errors in quadrature. 

Our best prediction is A th = (-3.0t^)10" 4 . It in- 
cludes B meson fragmentation effects (i.e. k — 0.84) and 
the values of the production fractions measured at the Z 
pole 3 [19(. For the sake of being more conservative, the 
uncertainty is derived by adding linearly the theoretical 
uncertainty quoted in Ref. [l[ with the uncertainty from 
B meson fragmentation quoted in table HI 



3 Note also that the world average value for fd is very close to the 
Please refer to sec. 3.1.3 of Ref. [Tj| for a detailed overview of one from the Z pole, while the one for f s is approximately the 

that difference and some of its implications. average of the two measurements figj] . 



IV. DISCUSSION 

It is obvious that the fragmentation effects analyzed 
in this paper cannot explain the large difference between 
the D0 measurement jjj and the SM prediction. More- 
over, additional sources of asymmetry are unlikely to be 
numerically significant (see appendix [A"|) . 

The most attractive possibility for explaining this ap- 
parent discrepancy is bSM physics. Clearly, the results of 
this paper will directly propagate in any analysis involv- 
ing new physics, as long as the effects of bSM physics are 
confined to the mixing of the B mesons, as expected. In 
the following we present the implication of the D0 mea- 
surement - including the fragmentation effects studied in 
this paper - for the extraction of ^4|l- 

On fig. Q] we present the D0 measurement in the ^4g L — 
Ag h plane. We give the existing constraints from in- 
dependent measurements on A$ L (from B-factories [l9j]) 
with a lig ht- g rey band, and on A gL (from the D0 mea- 
surement [20| of B® —> [i + D~ X) with a blue band. We 
plot the D0 dimuon asymmetry measurement in two 
ways. First we do not include fragmentation effects, i.e. 
K = 1, and take production fractions as measured at the 
Tevatron (with a dark-grey band); this is the same as 
figure 17 in Ref. [l|. Second, we plot the measurement 
including fragmentation effects, i.e. k = 0.84, and with 
production fractions as measured at the Z pole (with a 
red band). 

In all cases the plots are obtained by solving Eq. (|2"Tj) 
for Ag L and replacing with the measured value 
A exp = (-96 ± 29) 1CT 4 . In units of [1(T 2 ] we get: 

A* SL = aAl L +p. (22) 

The coefficients a and j3 are given in table HU The width 
of the two bands reflects the uncertainty in the measure- 
ment only. The end-lines of the bands on figure Q]can be 
obtained by adding ±0.796 (for the solid red band) and 
±0.588 (for the solid black band) to the coefficient f3 of 
the corresponding central line. 

TABLE II: The coefficients in Eq. ([22]). This paper is the 
center of the solid red band on figure [T]and the central red 
line on figure [2] D0 (same as figure 17 in Ref. y|) is the 
center of the solid black band on figure [T] and the black line 
on figure [2] Blue C,l,2 are the blue central and end-lines on 
figure Red 1,2 are the two red end-lines on figure [21 



Line 


This paper 


D0 


Blue C 


Blue 1 


Blue 2 


Red 1 


Red 2 


a 


-1.745 


-1.027 


-1.221 


-1.360 


-1.128 


-1.944 


-1.613 


P 


-2.636 


-1.946 


-2.132 


-2.266 


-2.043 


-2.826 


-2.509 



To get an insight into the sensitivity separately to frag- 
mentation (through k) and production fractions (through 
the values of f q ), on figure Owe plot the central values 
for the following scenarios: 
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1. (black line): k — 1 and production fractions 
A-(Tev), as in Ref. Q , 




-5 -4 -3-2-1012 



A d SL [lff 2 ] 

FIG. 1: The D0 measurement represented in the Ag L — ^IJl 
plane. The blue horizontal band is the constraint from the D0 
measurement [2g] of B° — >• fj, + Dj X; the light grey vertical 
band is the constraint from B-factories The solid black 
band is the one from Ref. [l| and the solid red band is our 
result. The small red circle is the SM value (uncertainty is 
not to scale). 




-5 -4 -3-2-1012 



A d SL [ia 2 ] 

FIG. 2: Variation due to fragmentation (i.e. k) and produc- 
tion fractions fi. Plotted are the three scenarios described 
in the text. Central values are shown. The black line corre- 
sponds to the center of the dark grey band on figure [1] and 
the central red line corresponds to the center of the red band 
on figure [T] 

The uncertainty on n corresponds to the ±5% uncer- 
tainty on A t h specified in table Q] We do not present the 
uncertainties from the measurement; they can be inferred 
from fig. [TJ 



6 



There are two important features that arise from our 
analysis. First, there is a point where the sensitivity 
to fragmentation functions/production fractions is zero. 
Strikingly, this point is precisely at the place where the 
three independent measurements agree with each other. 
An important implication is that the analysis in this re- 
gion is robust with respect to these parameters. 

Second, we note that each one of the two effects stud- 
ied in this paper: B fragmentation and sensitivity to 
the values of the B meson production fractions, pushes 
the band of the dimuon asymmetry measurement further 
away from the SM expectation in the region of small 
|-Ag L |. This result further strengthens the arguments 
given in Ref. [13] that new physics contributions are un- 
likely to only affect A %r , as one might naively expect; 
see also Refs. [28[ and [29(. Indeed, from figures ITI21 we 
conclude that the effects considered in this paper make 
it even more likely that the values of both A§ 
are significantly away from the SM expectations 



SL an d ^sl 



Note added. After this article was submitted for 
publication, a new measurement of the ratio f s /fd was 
performed by the LHCb collaboration [30]. The result of 
Ref. [30T ] is in complete agreement with the measurements 
°f fs.d at the Z pole. This new measurement confirms our 
QCD factorization-based arguments in favor of the val- 
ues of the fragmentation fractions as measured at the Z 
pole. Regarding figure [21 the measurement of Ref. [30l ] 
favors the set of red over blue lines. 
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Appendix A: Backgrounds 



are backgrounds and they only contribute to the denom- 
inator of Eq. j2j. In presence of such backgrounds the 
expression for the asymmetry can be cast in the form: 



_ ,(bkg=0) 



th 



1 + %uts) ' 



(5 (cuts) 



N ++ 



bk g 



Y++ , 

signal 



N 



signal 



N++ 
N++ , 

signal 



(Al) 



where j4[h ks °^ is the background-free asymmetry stud- 
ied in the previous sections. The correction S is estimated 
in Ref. to be ~ 0(10%). From Eq. (TAT]) it is evident 
that the effect of backgrounds is to decrease the value of 
A t h- The background correction is also cut dependent. 

Higher order effects can be an additional source of 
backgrounds. An example are partonic final states of 
the type bb + X or bb + X followed by two "right" decays. 
Such contributions are likely small, contributing at most 
few percent to A t h and likely much less. The reason their 
contribution is small is that such states are generated in 
processes suppressed by two powers of the strong cou- 
pling as, for example in the leading order contribution 
to the process pp —> bbbb. Alternatively, one can think of 
such processes as "usual" bb symmetric final states where, 
for example, the b quark perturbatively fragments into a 
6 .10] thus producing an inclusive final state with two b 
quarks. 

On the other side, the "gain" factor due to the absence 
of "wrong" type of decays is small for B s and Bd mesons. 
It can be estimated as follows. The suppression factor of 
a "wrong" decay with respect to a "right" decay is given 
by the ratio of time integrated probabilities (see |3.l3l|): 



a 



T[B q -> B q ] 



q T[B q -^B q ] 2 + x 2 ' 



(A2) 



Consider same-sign dimuon final states that contain 
or pairs with equal probability. Such events 



for q = s, d. Numerically C s « 1 and Cd ~ 0.2. 
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